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f^ ■ Abstract Using linear Vlasov theory of plasma waves and quasi-linear theory 

O^l I of resonant wave-particle interaction, the dispersion relations and the electro- 

^H ■ magnetic field fluctuations of fast and Alfvcn waves are studied for a low-beta 

multi-ion plasma in the inner corona. Their probable roles in heating and acceler- 
ating the solar wind via Landau and cyclotron resonances are quantified. In this 
paper, we assume that (1) low-frequency Alfven and fast waves, emanating from 
the solar surface, have the same spectral shape and the same amplitude of power 
spectral density (PSD); (2) these waves eventually reach ion cyclotron frequen- 
cies due to a turbulence cascade; (3) kinetic wave-particle interaction powers 
the solar wind. The existence of alpha particles in a dominant proton/electron 
plasma can trigger linear mode conversion between oblique fast- whistler and hy- 
<~i ' brid alpha-proton cyclotron waves. The fast-cyclotron waves undergo both alpha 

Mh. and proton cyclotron resonances. The alpha cyclotron resonance in fast-cyclotron 

Q ' waves is much stronger than that in Alfven-cyclotron waves. For alpha cyclotron 

f-j , resonance, an oblique fast-cyclotron wave has a larger left-handed electric field 

c/3 ' fluctuation, a smaller wave number, a larger local wave amplitude, and a greater 

, ^, , energization capability than a corresponding Alfven-cyclotron wave at the same 

wave propagation angle 9, particularly at 80° < 9 < 90°. When Alfven-cyclotron 
or fast-cyclotron waves are present, alpha particles are the chief energy recipient. 
1^ ' The transition of preferential energization from alpha particles to protons may 

be self-modulated by differential speed and temperature anisotropy of alpha 
particles via the self-consistently evolving wave-particle interaction. Therefore, 
^\i ] fast-cyclotron waves as a result of linear mode coupling is a potentially important 

^J • mechanism for preferential energization of minor ions in the main acceleration 

f-^ I region of the solar wind. 

Keywords: fast wave, linear mode coupling, solar wind 



1. Introduction 



C^ . The solar wind is a ubiquitous super- Alfvenic stream of magnetized charged par- 

ticles from the Sun, consisting of protons, electrons, and minor ions. According 
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to spectroscopic observations of extreme ultraviolet emission lines in the inner 
corona, solar wind ions exhibit a large kinetic temperature anisotropy, and minor 
ions are likely to be preferentially heated and accelerated ( |Kohl et ai, 1998 



Li et al., 1998[ [Cranmer et al., 1999). In the interplanetary space, minor ions 



have mass-proportional temperatures and similar bulk flow speed, whose drifting 
speeds relative to protons are at an order of local Alfven speed ( [Marsch et ai, 1982H 
[Neugebauer et ai, 1996 Goldstein et al., 2000 1 . And proton temperature in the 



direction parallel to the background interplanetary magnetic field is lower than 
that in the perpendicular direction. Double components of tenuous fast beam and 
dense core frequently occur in the velocity distribution function (VDF) of pro- 
tons, with a differential speed of 1.2 ~ 1.3 Alfven speed ( [Marsch et ai, 1982[ ). The 
combination of these remote-sensing and interplanetary in-situ measurements 
substantiates that some energy should be continuously deposited into the solar 
wind ions via some physical mechanisms, and thus power the hot coUisionless 
outflow from the Sun. 

The solar wind has been an outstanding puzzle for its heating and accelera- 
tion since the theoretical prediction ( [Parker, 1958[ ) and observation verification 
( [Neugebauer and Snyder, 1962P . A fraction of the mechanical energy in the Sun's 
internal convective motions must be converted into thermal energy above the 
photosphere. Ion cyclotron waves are generally considered as a promising mech- 
anism that can preferentially heat and accelerate minor ions. Ion cyclotron 
frequencies in the extended corona vary from 10 to 10^ Hz, whereas observed 
frequencies of pronounced oscillation on the solar surface are typically ~ 0.01 Hz. 
How strong enough ion cyclotron waves could be generated is still debated. In 
terms of the origins of ion cyclotron waves, there are two main viewpoints of base 
and secondary generations. For the base generation, high-frequency Alfven waves 
are directly launched during small-scale magnetic reconnection of micro-flares in 
rapidly evolving super-granular networks, then propagate up through the corona, 
gradually reach cyclotron frequencies of various ions, and are finally damped 
over a very short distance ( [Axford and McKenzie, 1992[ ). For the secondary 
generation, high-frequency ion cyclotron waves are gradually converted from low- 
frequency Alfven waves in the corona, and such conversion process is probably 
from a direct turbulence cascade ( Li and Habbal, 2001[[Markovskii et ai, 20"T0| , 



local nonlinear kinetic instabilities ( Markovskii and HoUweg, 2004[ ), and Debye 



scale electron phase space holes in electron Langmuir turbulence ( Cranmer and Ballegooijen, 2003P .| 



In this paper, we assume that (1) the plasma turbulence cascade is the energy 
source of the solar wind and (2) the solar wind heating is due to Landau and 
cyclotron resonances with a continuous turbulence spectrum. 

The solar wind turbulence is an ensemble of fluctuations with random phases 
and a broad range of wave vectors. Such fluctuations arc manifested in mag- 
netic field ^B, electric field (5E, plasma bulk flow velocity SV, and so on. The 
fluctuation energy is usually transferred through a forward cascade. Namely, 
the energy is injected from a preexisting population of magnetohydrodynamic 
(MHD) waves at very low frequencies / and wave numbers k, nonlinearly trans- 
ported through successively shorter wavelengths via wave- wave interaction, and 
eventually dissipated at ion kinetic scales via coUisionless wave-particle interac- 
tion. The evolution of 5B and 6'V in the presence of a background magnetic 
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field Bq can be represented in the Elsasser variables Z^ = SY ± 6'B/^/A^Fpo as 
follows: 

i^T(^A-V)Z± + (ZT.V)Z±=-VP + i^, (1) 



k 



Bo > k^bx 
Bq w k±bx , 



weak turbulence, 
strong turbulence. 



(2) 



Here va — Ba/^'inpo is the Alfven speed, po is the mass density, P is the 
pressure that is determined from the incompressibility condition of V • Z* = 0, 
.F is a large-scale forcing, fcy and k±_ are parallel and perpendicular wave num- 
bers, and b\ (^ Bq) is oscillation amplitude at the scale A « l/fc_L- Z* is 
imbalanced at a macro scale, as supported by nonzero cross-helicity in terms of 
SV ■ (5B (Goldreich and Sridhar, 1995 He et ai, 2011). There is more power in 



Alfvenic fluctuations travelling away from the Sun than towards it. Even if bal- 
anced overall, MHD turbulence is always locally imbalanced in creating patches 
of positive and negative cross-helicitics. The conservation of cross-helicity re- 
sults in a hierarchical structure of MHD turbulence. In the corona, turbulence 
amplitude is constrained by interplanetary scintillation (IPS) observations of 
electron density spectra ( |Cole, 1978 ). Interplanetary in situ measurements of 
turbulence (5B as a function of /, fcy, and k± have shown a power law spec- 
trum of magnetic field fluctuation, and a large anisotropy in favor of quasi- 



perpendicular wave propagation (Matthaeus et al., 1990 Leamon et ai, 1998 



[Dasso et ai, 2005|[Ga"ry and Smith, 2009| pacBride et ai, 2010 1 



\5B{f)\'^.r 



(3) 



|<5B(A.||)P = ^|<5B(k)p«fcp 

\sBik^)f = j2\^Bik)f^kr' 



(4) 

(5) 



ail > a I 



(6) 



The solar wind turbulence in an inertial range of 10"''' < / < 0.2 Hz is similar to 
the classical Kolmogorov picture of fluid turbulence, and the spectral index a in 
Equation ([3]) is approximately 5/3 ( [Leamon et ai, 1998[ [Horbury et ai, 2005 ). 
Meanwhile, the power spectral density is anisotropic with respect to Bq, as indi- 
cated from ay ss 2 in Equation (g]) and a± w 1.6 in Equation ^ ( [Gary and Smith, 2009[ ).| 
Moreover, within the inertial range, the Elsasser variables Z^ were also found 
to be anisotropic and have two distinct scaling subranges ( |Wick et al, 20ll| ) . A 
transition between weak and strong turbulence regimes (Equation ([2])) is feasible 
by changing a large-scale forcing F (Equation llj), as demonstrated from numeri- 
cal simulation of anisotropic incompressible MHD turbulence ( [Perez and Boldyrev, 2008] ) .[ 
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The nature of the MHD turbulence cascade Hes in the interaction of wave pack- 
ets moving with Alfven velocities. In the framework of the weak turbulence 
scenario, interacting wave packets are very large such that the original wave 
packets decay before they pass through each other. Because the auto-correlation 
and cascade time scales are always of the same order of magnitude in the 
weak turbulence theory, any model that implies a large number of collisions 
among wave packets for an efficient energy cascade should belong to the strong 
turbulence theory ( [Gogoberidze et ai, 2009[ ). In the strong turbulence, interac- 
tion is strong, resonance conditions are not in force, hence a parallel cascade 



is possible ( Perez and Boldyrev, 2008 ). The forward cascade of (5B could be 
modelled as a combination of advection and diffusion in the wave number k 
space dCranmer and Ballegooijen, 2003D . Equations (JH) and (O for turbulence 
anisotropy in the inertial range could be derived by a conjecture of dynamical 
alignment that polarizations of i5B and JV tend to align during the turbulence 
cascade towards successively smaller scales ( [Goldreich and Sridhar, 1995] ). A 
good fraction of the power spectral density (PSD) observations is well interpreted 
by the current turbulence theories. 

The turbulence in the dissipation range is poorly understood. Spectral prop- 
erties in the dissipation range almost certainly depend upon the details of the 
damping, though those in the inertial range are relatively robust and independent 
of detailed turbulence cascade processes. The observed wave number at a break- 
point, which separates the inertial and dissipation ranges, may scale as inverse 



proton inertial length kc/ujpp ~ 1 (Leamon et al., 1998 Smith et at, 2001) or 



cyclotron radius kvp/flp ~ 1 ( [Bale et al, 2005 1, with Wpp, Op, and Vp denoting 



proton plasma frequency, proton cyclotron frequency, and proton thermal speed 
respectively. Above the inertial range, the PSD has a distinct breakpoint at 
0.2 < / < 0.5 Hz; At higher frequencies of / > 0.5 Hz, the PSD is steeper 
with a broader range of 2 < a < 4 in Equation ^ ( [Leamon et al, 1998} 



[Gary and Smith, 2009[ ). For instance, ISahraoui et'oU ([MJ^)) found that a PSD 
above the inertial range consists of two distinct regimes with successively larger 
values of a in Equation (|3]); a is about 2.5 from 0.4 to 35 Hz for a dispersion 
cascade and 4 from 35 to 100 Hz due to strong wave damping. Furthermore, it 
is still unclear what are the constituent wave modes of the solar wind turbu- 
lence. The macro-scale fluctuations in the solar wind and corona are Alfvenic in 
nature ( [Barnes and HoUweg, 19741 ). In the MHD regime, Alfven waves can be 
converted into fast waves by density fluctuations, and vice versa through resonant 



three- wave interaction (Chandran, 2008). At the kinetic scale, highly oblique 



whistler waves exist according to the observed properties of large magnetic 
compressibility in the solar wind ( [Gary and Smith, 2009[ ) . Kinetic Alfven waves 
are identified from the observations of largely enhanced electric fluctuation spec- 
trum and super- Alfvenic wave phase speed at scales of ion thermal gyro-radius 
( [Leamon et ai, 1998[ [Bale et al., 2005[ ). Parallel propagating Alfven-cyclotron 
waves are possibly present as evidenced by the measured angle distribution of 
magnetic helicity ( [He et ai, 201 ip . These observations have revealed the impor- 
tance of quasi-perpendicularly propagating waves in the heating of the solar 
wind. 

Increasing attention has been paid to oblique propagation owing to the nat- 
ural occurrence of oblique waves from refractions in the inhomogencous solar 
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wind, production by density gradients, and generation from MHD turbulence 

dLi and Habbal, 2001[|Gary and Nishimura, 2004{|Osmane et al, 2010[[Li et ai, 20T0\ § 



|Markovskii et ai, 2010 ). Indeed, for the Alfven-cyclotron wave branch, paral- 
lel propagation has its inherent difficulty to reconcile both the heating and 
the acceleration of the solar wind ([Hollweg and Isenberg, 2002). The nature 



of obliquely propagating waves is subject to debates ( Gary and Smith, 2009[ 



Podesta et ai, 2010 1. According to interplanetary measurements, these constituent! 



waves may be kinetic Alfven waves ( Bale et ai, 2005[ [Sahraoui et al., 2010[ ) or 



fast- whistler waves ( [Gary and Smith, 2009 1. In a warm plasma such as the in- 
terplanetary solar wind, oblique Alfven waves can generate a fast proton beam 
through Landau resonance ( |Li et ai, 2010[ |Osmane et ai, 2010[ ) . In a cold pro- 
ton/electron plasma, quasi-perpendicular fast-whistler waves break up into Bern- 
stein modes near the first few harmonics of proton cyclotron frequency ( |Li and Habbal, 2001H 
|Markovskii et ai, 2010[ ). Due to the existence of minor ions in a solar corona- 
like plasma, oblique fast-whistler waves are linearly coupled with hybrid cy- 



clotron waves ( Li and Habbal, 2001 ) . Depending on the abundance and relative 
speed with respect to protons, minor ions may change wave dispersion rela- 
tion dramatically, induce a transition of wave polarization from right- to left- 
handed sense, and can absorb more energy from oblique fast-cyclotron waves 
than Alfven-cyclotron waves for intermediate angles 30° - 50° of wave propaga- 
tion ( |Li and Habbal, 2001[ ). However, Alfven-cyclotron and fast-cyclotron waves 
in the inner corona have few empirical constraints for their exact properties 
such as generation mechanisms, propagation modes, and power levels. With 
favorable plasma parameters, both Alfven and fast waves can be potential can- 
didates for preferential energization of minor ions in the solar wind. On the 
basis of current understanding, both fast-whistler and Alfven-cyclotron waves 
at both quasi-parallel and quasi-perpendicular propagation angles probably con- 
tribute to kinetic-scale turbulence in the solar wind ( [Gary and Smith, 2009 
|He et al., 20TT1 ). 

In this paper, both Alfven and fast waves are investigated to explore their 
roles in energizing the low-beta multi-ion solar wind plasma in the inner corona. 
As a continuation to the work of lLi and Habbal1(|200ip . the linear mode conver- 
sion between fast-whistler and hybrid cyclotron waves is further quantified by a 
detailed parametric study of various plasma parameters. A quasi-linear theory 
of resonant wave-particle interaction ( [Marsch and Tu, 2001[ ) is used, and we are 
able to study the parallel and perpendicular heating of solar wind species due to 
Landau and cyclotron resonances. We give our model in Section [2j and describe 
wave dispersion relations and accompanying electromagnetic field oscillations in 
Section [21 Then the roles of wave propagation angle, minor ion presence, relative 
flow of minor ions, and plasma temperature are successively analyzed in Sections 
miilini and [71 We discuss the ensuing nonlinear stage of wave-particle interaction 
in Section [SI and summarize our results in Section [HI 



2. Methods 

To evaluate the energization effect of kinetic wave-particle interaction on the 
multi-species solar wind, two methods are sequentially combined in this pa- 
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per. First, using linear Vlasov theory ( |Stix, 1992} |Li and Habbal, 2001[ ), the 
properties of Alfvcn and fast waves such as their dispersion relations and elec- 
tromagnetic polarizations are studied. Second, using a quasi-linear theory of 
resonant wave-particle interaction ( [Marsch and Tu, 2001[ ), the previously de- 
rived wave properties arc adopted to calculate the acceleration and heating 
rates of different plasma species. With the two-step approach, the wave-particle 
interaction can be quantified to demonstrate the most viable and efficient way 
for the solar wind energization. 

Since we are interested in the main solar wind acceleration region, typical 
plasma conditions in the inner corona are chosen in our analyses. Given a 
background magnetic field Bq, a wave vector k is decomposed into a parallel 
component fcn and a perpendicular one k±. The propagation angle of a wave 
= arctan(fc_L/A:||) is defined with respect to Bq. Our Cartesian system of 
(x, y, z) coordinates is constructed with Bq along the z-axis and k within the 
xz-plane. The symbols and corresponding denotations used in this paper are 
summarized in Table [T] In this paper, only electrons c, protons p, and alpha 
particles a are considered in the solar wind. For simplicity, VDF of each species 
is approximated to be bi-Maxwellian. The parallel (perpendicular) thermal speed 
vuj {v±j) and gyrotropic VDF fj(vu,vj_) are expressed as follows: 



/2fcBT||,/m, , (7) 

v±j = y2kBT±j/mj , (8) 

/,(v) = /,(«,,, «^) = ^ , exp f- ^"ll">-^' - 4-] ■ (9) 

The dispersion relation ui = w(k) and its inherent electromagnetic field {SE, 
i5B) of a plasma wave are calculated using linear Vlasov theory ( |Stix, 1992{ 
|Li and Habbal, 2001[ ). The complex wave frequency a;(k) = Wr(k) -I- iujiik), 
(5B(k), and (5E(k) are found as a function of wave vector k. The damping rate 
|wi| generally increases with wave number fc. |Gary|p993p suggested that heavily 
attenuated waves are described by a rough approximation: 

c^i(k) < -^ . (10) 

Z IT 

and that only weakly damped waves can persist. Given (5E(k), polarization P 
and unimodular electric vector e£;(k) are given by 

P = *T^^ , 11 

(5Ey IWrl 
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Consequently, the magnetic field fluctuation i5B(k) and its power spectrum 
|(5B(k)|2 are 

<5B(k) = ^kx<5E(k), (13) 

\SB(k)f = (^^) (l-|e,-eMk)P)|«-E(k)p, (14) 

The wave-particle interaction via Landau and cyclotron resonances is efhcient 
in transferring wave energy to particles in a coUisionless plasma. With the in- 
troduction of Doppler-shiftcd frequency w'(k) = a;i.(k) — fcii t/iij, the nth-order 
resonant parallel speed in the bulk frame of species j is defined as: 

II I 

Then, the resonant factor ^j for species j is defined as: 

^j{K'n) = — = ^ , n = 0,±l,±2,... . (17) 

Here n = corresponds to the Landau resonance, n ^ the cyclotron resonance. 
Only the Landau resonance ^° and the first-order cyclotron resonance ^~'" may 
appear for the plasma parameters prescribed in Tables [5] and [31 

The acceleration and heating rates of a species due to any wave mode were 
formulated bv lMarsch and T"ul(|200ip using a quasi-linear theory, and arc applied 
in this paper. In the quasi-linear theory, different waves are simply considered 
to be linearly superposed without any interference. Specifically, for species j, 
its parallel acceleration rate ay, parallel heating rate Q\\j, and perpendicular 
heating rate Q±j are integrated over various wave modes M and all resonance 
orders n, and are listed below: 



%/-^ - E Ur h^ T-^-^7n:V^^^(k,n) 2fc||..,(k,n) 



Q^,/m,) «f"ooV2^/ \k) l-|ek-e£(k)|- ^ ^^^ 

(19) 

Herein, the resonance function Rj (k, n) gives energy absorption of species j 

from a wave mode |5B(k)p. An electric field oscillation could be factorized into 

a left-handed polarization and a right-handed one. The factorized components 
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of unimodular electric field e£;(k) and total electric field (5E(k) are expressed as 

e|(k) - e£,(k)±ieBy(k) , (20) 

6E^{k) = 5E^{k) ± i 6Ey{k) . (21) 

In Equations (POJ) and (PT|) . subscripts + and — denote left- and right-handed 
polarizations, respectively. The transverse wave number k± is compared with 
the ion cyclotron radius of thermal motions v±j /flj , and hence described by a 
parameter Xj. 



2 V ^j 
Further, the resonance function i?j(k, n) is defined as 

Rj{k,n) ^ 2Nj{k,n) cxp{-Xj)[a^I,,^i{Xj)+a+In+i{Xj)+aoIn{Xj)] , (23) 

N,{k,n) = 0FiexpK-(k,n)]x(o(k,n)|^+n^:^| , (24) 



|2 



II L ^j\\ «'ll'^lb- 

ao{k,n) - 2 ^-f ^ \eE,\^-2Im{e},yeEz}^^^ '- + 2 X, \eEy 

(25) 
a^(k,n) = ±^ |e±p ± Rc{e^.e^}^ "('^^^'^ ""^ T A, Im{e^,e|} . (26) 

Here / is a modified Bessel function. For parallel propagation, the above form 
of Rj (k, n) can be largely simplified as 



i?j(0,fc||,ri) =iVj(0,fc||,n) 



,+ |2^A .l„-|2 ... ,„ ,2/'^(fc||) 



,1 |ej;r + ^n,-i le^r + 4(5„,o |e£;^ 



•|| «ij . 

(27) 
Note that the resonance function Rj (k, n) essentially depends on the resonance 
factor ^j and the square of the re-scaled perpendicular wave number Xj = 
{k^vp^/^j)"^ /2 as exponential functions of exp(— ^|) in Njik^n) (Equation ([24])) 
and exp(— Aj) in Rj{k,n) (Equation (051) ') respectively. As ^j or Xj is increased, 
exp(— ^^) or exp(— Aj) dramatically declines to zero, and energization effect 
quickly becomes negligible. 

We assume that plasma turbulence |(5B(k)p is the ultimate energy source 
of the solar wind. For simplicity, energization contributions of Alfven and fast 
waves are individually evaluated in terms of the energization rates an, Qii, and 
Q± in Equation ([TO)) . Both Alfven and fast waves are assumed to have the 
same PSD |(5B(k)p at the same wave vector k. A specific form for |(5B(k)p in 
Equation (J18p is beyond the scope of this paper. Initial prescriptions of various 
plasma parameters are given in Tables [2] and [3] Then the values of oii, Qn, and 
Q±_ (Equation ([TO)) ) normalized to |5B(k)p (Equation ([T5)) ) are derived by a 
combined two-step approach ( |Li and Habbal, 2001[|Marsch and Tu, 2001[ ). Here 
a detailed parametric study is conducted to analyze how cascaded turbulence 
energy is repartitioned among multiple constituent plasma species. 
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3. Dispersion Relation and Electromagnetic Field Oscillation 

Dispersions of Alfven and fast waves in a proton/alpha/electron plasma are 
illustrated for various propagation angles 6 in Figures[T][21 and|31 In the presence 
of alpha particles, an Alfven-cyclotron wave branch of < Wr < Vt-p is generally 
split into two parts near the cyclotron frequency of alpha particles Vta- The lower 
branch of < a;,. < Via corresponds to Alfvcnic alpha cyclotron waves, and the 
upper one of Via < Wr < fip refers to hybrid alpha-proton cyclotron waves. 
Obviously, the alpha cyclotron resonances occur at much smaller wave numbers 
for the hybrid cyclotron waves than the Alfven-cyclotron waves. Though having 
not been measured experimentally, kinetic hybrid waves with uj-,- « Via and fc « 
may be locally generated in a cold multi-ion plasma as a result of linear mode 
coupling with fast-whistler waves ( |Li and Habbal, 2001[ ). Avoiding a crossing 
between the uj — k dispersion curves of hybrid ion-ion cyclotron and fast-whistler 
waves, fast-whistler waves are possibly transformed into fast-cyclotron waves. By 
comparing Figurc[T^ with[T]:, one can see that the dispersion of oblique fast waves 
is dramatically changed by the presence of alpha particles in a dominant elec- 
tron/proton plasma. Fast-whistler and alpha-proton cyclotron waves exchange 
their identities at alpha cyclotron frequency ui-c = JIq, — 0.5 ilp. Such linear 
mode conversion in Figure [TJ; occurs when fast waves are oblique enough, i.e. 
the propagation angle 9 > 31° in this paper (Figure [5^). Once mode conversion 
happens, the conversion point in the dispersion curve is basically at one identical 
point of w,./ rip = 0.5 and k v^/ ftp = 0.5, irrespective of wave direction 9. At a 
sufficiently long wavelength, fluctuations of both Alfven and fast waves are essen- 
tially undamped. For ion cyclotron waves, as the parallel wave number fcii reaches 
a characteristic dissipation value k^ , relevant ions are in cyclotron resonance and 
the damping begins abruptly. In a proton/electron plasma of 0.001 < /3p < 0.1, 
kd ~ flp/vA for the proton cyclotron resonance was given by |Gary and Borovsky| 
(|2004p : In a cold alpha/proton/electron plasma, fed ~ ^a/vA = 0.5rip/z;A for 
the alpha cyclotron resonance is found in this paper. Specifically, for nearly 
parallel Alfven waves oi 9 = 1°, the damping rate |wi/r2p| is increased from 
2 X lO"'^ at fc ua/ rip = 0.45 to 4x10^3 at fcwA/f^p = 0.6 (Figurefl})) . Meanwhile, 
damping of fast-whistler waves is only pronounced at harmonic ion cyclotron 
frequencies for oblique propagation ( |Li and Habbal, 2001 1. However, due to lin- 



ear mode conversion, fast-cyclotron waves will experience a dramatic increase 
of damping rate |a;i|. as shown in Figure [TJi. Heavily damped waves disappear 
swiftly, as described in Equation (|10p ( [Gary, 1993[ ). Approaching ion cyclotron 
resonance frequency, smaller parallel wave numbers fcy and larger total wave 
numbers k are found as 9 increases to be more oblique, especially at w 90° 
(Figures [2] and |3]) . When a turbulence cascade is present, Alfven-cyclotron, fast- 
whistler, and fast-cyclotron waves may continuously transfer wave energy from 
very low frequencies along their respective dispersion curves, and ultimately 
become constituent modes of kinetic turbulence at ion cyclotron frequencies. 

Electromagnetic behaviors of Alfven and fast waves are displayed in Figures 
EJ [5l and |6l As fast-cyclotron and Alfven-cyclotron waves become more oblique, 
phase speeds uj,-/k involving cyclotron resonances are gradually reduced towards 
0, and electric field oscillations (5E become more aligned with wave vector k. 



SOLA: SoPh2011Rev2.tex; 25 February 2013; 23:47; p. 9 



Solar Physics 



The values of uj^/k and (1 — |efe • esj) determine the ratio of electric to magnetic 
field fluctuations |(5Ep/ |(5Bp (Equation dH])). For Alfven-cyclotron waves with 
Q° <e < 90° (Figures SJa, and[5]a) and fast-cyclotron waves with 31° < 6* < 90° 
(Figures HJ; andlHh.), |<5Ep/|(5Bp increases when becomes larger. At oblique 
enough angles, both fast and Alfvcn waves are nearly electrostatic and linearly 
polarized as a result of (5E being almost parallel to k. In particular, nearly- 
perpendicular fast waves become almost electrostatic, as |5Ep/(|(5Bp u^) could 
be as high as 50 at fcwA/ i^p = 3 (Figures |3^ and|B^). As a result, electric field 
fluctuations |(SE| of fast waves are significantly enhanced over their corresponding 
magnetic field fiuctuations |(5B|. As for Alfven waves, |(5E| and |(5B| • wa are 
equally strong for quasi-perpendicular propagation (Figures H^ and[3i), though 
|(5E| is far less than |(5B| • ua at large k for parallel propagation (Figures |3J 
and [5^). Meanwhile, for oblique propagation, parallel electric field fiuctuations 
\5E\\\ are enhanced over the whole spectrum. The overall sense of the electric 
field polarization is left-handed for parallel Alfven waves (Figures |3] and[5}3), 
right-handed for parallel fast waves (Figure [Sb): and nearly linearly polarized 
for perpendicular Alfven and fast waves (Figures IDd, Hf, [S}d, and EJd). Obvi- 
ously, linear mode conversion from fast waves to alpha cyclotron waves near 
alpha cyclotron frequency 51q abruptly increases the left-handed electric field 
oscillation (Figure ^p) . When present in the inner corona, both Alfven and fast 
waves could heat plasma species through Landau and cyclotron resonances. If 
a resonant factor \^j \ is large, only particles in the remote tail of the VDF of 
species j are in resonance with the wave. In linear Vlasov theory, a resonance 
factor of l^jl < 3 is a necessary but not a sufficient condition for non- negligible 
wave-particle interaction ( |Gary and Borovsky, 2004] ). Both |^°| < 3 and parallel 
electric field oscillation \5E\\\ > are required for noticeable Landau resonance; 
Similarly, both |^^| < 3 and left-handed electric field oscillation \SE^\ > 
correspond to cyclotron resonance. Due to nearly zero \SE\\ \ accompanying quasi- 
parallel waves. Landau resonances are very weak even with the synchronous 
condition of wave and particle phases, i.e. |^°| < 3. Conspicuous streaks of ^+ 
from 0° < 6* < 31° (FigurelBJ;) and ^+ from 0° < 6* < 90° (FigureHf) correspond 
to locally enhanced damping of fast waves at ion cyclotron frequencies ila and 
Up. However, as explained by ILi and Habball (j200ip . the localized cyclotron 
damping of fast waves by species j is unlikely to be important because of its 
localized occurrence in a very narrow frequency band across flj . Through Landau 
and/or cyclotron resonances, different plasma species absorbs energy from the 
turbulence power spectrum |(5B(k)p accordingly. 



4. Effect of Wave Propagation Angle 

The energization of any plasma species j due to turbulence PSD |(5B(k)p is 
estimated by quasi-linear transport coefficients (aiu-, Qiij, and Q±j) in Equation 
P^ and visualized in Figures [7] and |S1 Alfven and fast waves at each k are 
assumed to be equally strong; the magnitude of fluctuation |^B| is the same 
for both waves. The parallel acceleration a\\j, parallel heating Q||j, and per- 
pendicular heating Qj^j are basically determined by the resonance function R 
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in Equation (P5)) . Different species show different behaviors in terms of their 
energization from wave damping. Affven-cyclotron waves essentiahy do not in- 
teract with protons, because their asymptotic frequency fla is far away from the 
proton cyclotron frequency J7p (Figure[2^). Fast-whistler waves with propagation 
angles 0° < 6' < 31° are essentially undamped, as shown in Figure [8] Notice- 
able energy recipients of Alfven waves are alpha particles when 0° < 9 < 80° 
and electrons when 80° < 9 < 90° (Figure [7]); those of fast waves are alpha 
particles with 31° < 6* < 90°, protons with 31° < 6* < 80°, and electrons with 
80° < 9 < 90° (Figure [H])- The alpha cyclotron waves in Figure [7] have similar 
energization effects as the proton cyclotron waves in Figure [3 though being 
damped more heavily. The strongest energy recipient among all plasma species is 
the alpha particles via alpha cyclotron resonances with fast-cyclotron waves when 
31° < 9 < 90° (Figures [5^-c). Such an alpha energization from fast-cyclotron 
waves is increasingly significant, as the waves become more oblique towards 
9 = 90°. Moreover, for effective energization from fast-cyclotron waves, alpha 
particles occupy a lower part of the k\\ — 9 parameter space than cither protons 
or electrons. Based on the turbulence cascade scenario, the PSD is absorbed 
by resonant alpha particles from fast-cyclotron waves at lower frequencies and 
wave numbers, and hence has larger fluctuation amplitudes |5Bp. As a result, 
the combination of the largest |(5Bp in Equation ([T5|) and the largest (an^, Q||q,, 
QA_a) in Equation (J19p are responsible for the strongest alpha energization. In 
addition, compared with alpha cyclotron frequency fio,, frequencies uj^ at the 
alpha cyclotron resonance arc slightly higher for fast-cyclotron waves (Figure [TJ;) 
and slightly lower for Alfven-cyclotron waves (Figures [T^ and[T]:). Consequently, 
the first-order parallel resonant speeds WQ,(k, 1) are positive for fast-cyclotron 
waves and negative for Alfven-cyclotron waves (Equation (J16p ). According to 
Equations ([T5| and (fTO)) . fast-cyclotron and Alfven-cyclotron waves are respon- 
sible for parallel heating and cooling of alpha particles respectively. With Alfven 
waves, it is much easier for the VDF of alpha particles to develop into a kinetic 
anisotropy of T±a > Th^, and then trigger an ion cyclotron instability of alpha 
particles. In a cold plasma of /3p = 0.01, alpha particles arc the primary energy 
recipient from the damping of both Alfven and fast waves. This preferential 
energization of alpha particles is much more significant by oblique fast-cyclotron 
waves via linear mode conversion from their corresponding fast-whistler waves. 



5. Influence of Minor Ion Presence 

Without alpha particles in a proton/electron plasma, oblique fast- whistler waves 
can still be linearly converted to Bernstein waves at harmonic proton cyclotron 
frequencies of Qp and 2flp ( |Li and Habbal, 200 1| [Markovsku e^ ai, 20T0l ). At 
a perpendicular angle oi 9 = 90°, Bernstein waves are purely electrostatic, 
do not damp, and cannot affect ions ( |Stix, 1992] ). When 9 deviates from 90°, 
Bernstein waves are strongly damped. With linear mode coupling, highly oblique 
fast-Bernstein waves are strongly dissipated by protons, as previously demon- 
strated from linear Vlasov analyses ( |Li and Habbal, 2001[ ) and hybrid simulation 
dMarkovskii et al., 2010[ ) of fast-wave cascade spectrum. With the absence of 
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alpha particles, energization rates of protons and electrons are estimated by 
Equation ([TO)) . Only differing in the abundance of alpha particles n^ = 0.02 Tip 
(Tables [5] and [3]) , the two sets of initial plasma parameters in Sections [3] and 
[5] result in dramatically different energizations. Linear mode conversion from 
fast to Bernstein waves occurs, when the fast waves are oblique enough, i.e., 
propagation angles 6 > 64° . Specifically, such a conversion point of fast-Bernstein 
wave branch is at 2 ftp for 6A° < < 87° and Qp for 87° < 6* < 90°. Energization 
of protons and electrons arc noticeable for 64° < 6* < 85° and 80° < 6* < 90° , 
respectively. Because of the linear mode coupling, fast-whistler waves arc in- 
tercepted at frequencies Wj ~ fla for alpha cyclotron waves and uji ~ flp (or 
2 rip) for Bernstein waves. Parallel electric field \SE\^\ between fl^ and fip of 
fast-Bernstein waves are responsible for a much stronger electron energization 
in contrast with fast-cyclotron waves. Protons undergo larger energization by 
fast-Bernstein waves than by fast-cyclotron waves, particularly in perpendicular 
heating Q±p. However, fast-cyclotron waves primarily energize alpha particles, 
and such a strong ion energization cannot be found in fast-Bernstein waves. 
When alpha particles are present, fast waves are coupled to cyclotron waves 
instead of Bernstein waves as a result of the turbulence cascade. Approaching 
alpha cyclotron frequency fla (instead of fip), a fast-cyclotron wave has a smaller 
frequency Wr, a smaller wave number k, and larger wave energy density |(5B(k)p 
based on the turbulence cascade scenario. Moreover, the angle threshold for 
linear mode coupling is low for fast-cyclotron waves {9 > 31°) and high for fast- 
Bcrnstein waves {6 > 64°). Therefore, the existence of alpha particles could lead 
to more efficient damping of fast waves. 



6. Role of Relative Flow of Minor Ions 

Wave-particle interaction in a multi-ion plasma may be significantly modulated 
by a differential speed U^^j of minor ion j. Here U^^j is defined as the bulk flow 
speed relative to protons. The dependence of energization on C/||q, in a pro- 
ton/alpha/electron plasma is presented in Figure IH] for Alfven waves and Figure 
[Tn]for fast waves along propagation angle 9 = 60°. Initially cold alpha particles 
are swiftly accelerated and heated by both Alfvcn-cyclotron and fast-cyclotron 
waves. Fast-cyclotron waves resonate with alpha particles at frequencies Wj ~ fla 
and protons at uji sa f2p, and energize alpha particles at a low wave number 
k (Figures [TOk-c') and protons at a high k (Figures ITOH-f). Obviously, alpha 
particles are energized much more strongly than protons for Alfven-cyclotron 
waves at C/y^ < 0.2 wa (Figure IH]) and for fast-cyclotron waves at C/||q < 0.6 wa 
(Figure[TU]). When C/||c reaches 0.2 va, Alfven-cyclotron waves at 9 = 60° begin 
to resonate with protons, and alpha particles lose the initial energization advan- 
tage. Meanwhile, oblique fast-cyclotron waves continue to preferentially energize 
alpha particles. The appearance of fast-cyclotron waves occurs for 6 > 31° at 
Uiia/vA = (Figure 0) and 9 > 60° at U\ia/ va = 0.6 (Figure [lOj) . Obviously, 
as a response to continuously increasing C/iu, the minimum threshold of wave 
propagation angle 9 at which the coupling between fast-whistler and hybrid 
alpha-proton cyclotron waves occurs is increased accordingly. As a result of 
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continuous preferential energization, alpha particles gradually accumulate high 
enough t/iu and temperature anisotropy T^a/Tn^. When alpha particles lose 
resonance, the major wave-particle interaction is proton cyclotron resonance for 
Alfven-cyclotron waves (Figures |HJi-f ) and electron Landau resonance for fast- 
whistler waves (Figure fTOh). Productions of much higher C/||q and Tiq/T||q 
are inhibited, as excessive free energy of alpha beam and kinetic anisotropy in 
the VDF of alpha particles would be released to excite Alfven-cyclotron waves 
(Figures [5^-c). The self-consistent evolution of alpha particles in the parametric 
space of C/|u and Tlq/Thq could be further investigated by particle simulations. 
Here from Vlasov analyses of various plasma parameters, both fast and Alfvcn 
waves arc likely to contribute to the eventually sustainable differential speed and 
temperature anisotropy of minor ions in the solar wind. 



7. Dependence on Plasma Temperature 

Fast-cyclotron waves only exist in a cold plasma, because linear coupling of 
fast-whistlcr and hybrid-cyclotron waves disappears for large plasma beta. Para- 
metric studies of proton beta /3p are performed for Alfvcn and fast waves along 
propagation angle 9 = 60° (Tables [5] and [3]). Generally, in an increasingly 
hot plasma, waves are increasingly damped, and their effective uj^ — k disper- 
sion curves are consequently shortened according to Equation (jlOp . Max.(fc||) 
generally decreases with an increasing /3p, which is estimated at each f3p in 
the ku — /3p parameter domain according to weak damping condition (Equa- 
tion (|10p ). Exceptional increases of Max.(fc||) are found from a transition of 
— 1.2 < logjo(/3p) < —1 and a jump across logig(/3p) = —0.95. Such exceptions 
in the trend of Max.(/s||) versus /3p are due to a significant dispersion change 
of corresponding wave branch. As /3p increases, the dispersion of Alfven waves 
finally changes to cross the alpha cyclotron frequency fia, and fast waves ulti- 
mately cease the conversion to hybrid alpha-proton waves. The sudden change 
of wave dispersion corresponds to a certain value of /3p. The critical point is 
logio(/^p) = —1-2 for Alfvcn waves and log]^o(/3p) = —0.95 for fast waves. Across 
the critical points, wave branches extend to higher frequencies Wr, and survive 
for higher wave numbers fcy . Alpha particles are the chief energy recipient from 
both Alfven-cyclotron waves and fast-cyclotron waves (logio(/3p) < —0.95). The 
energization effects of Alfven-cyclotron waves weaken as /3p increases, while those 
of fast-cyclotron waves peak in the range of —2.2 < log]^o(/3p) < —1.2. In addi- 
tion, without mode conversion for log2o(/?p) — ^0.95, oblique fast-whistler waves 
not only energize ions around ion cyclotron frequencies {fla and Hp), but also 
dissipate a noticeable portion of wave energy to electrons in the form of parallel 
heating <5||o- To the energy absorption of fast- whistler waves, the contribution of 
electrons is comparable to that of protons and alpha particles in a warm plasma. 
Only in a cold plasma, linear coupling of fast and cyclotron waves exists, and the 
electron heating could be negligible for damping of fast-cyclotron waves. Hence, 
in the inner corona with /3p f* 0.01, oblique fast-cyclotron waves are expected 
to play a more important role in energizing minor ions, provided that they are 
equally as strong as Alfven-cyclotron waves. 
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8. Discussion on Nonlinear Stage of Wave- Particle Interaction 

It is still an enigma how the macro-scale heating and acceleration of the solar 
wind are powered by waves on ion kinetic scales and how these waves are gener- 
ated by multi-scale coupling. The solar wind turbulence becomes fully developed 
on a time scale that is very short, compared to the expansion time of the solar 
wind along open magnetic flux tubes. So the turbulence evolution in the wave 
vector space could be considered to develop in a spatially homogeneous plasma, 
so we can adopt spatially uniform energization rates (Equation ([T5)) ') in this paper 
(see also Cranmer and Ballegooijen (j2003|) ). In Equation ((19)) . a\\j is probably 



not an important contributor to the solar wind acceleration, as the magnetic 
mirror force due to anisotropic heating of Q\\j and Q^j is typically stronger. A 
plasma species j is energized by dissipation of its ion cyclotron waves. In a frame 
moving with the phase speed of the waves, a diffusion in the VDF of species j 
(Equation ([9])) occurs along contours of constant energy. Such a diffusion process 
is faster for species j with a lower charge-to-mass ratio. This wave-particle 
interaction was formulated into a kinetic shell model Psenberg, 2001[ ), which 
quantitatively describes ion energization by bidirectionally propagating ion cy- 
clotron waves. As a result, the VDF in Equation ([9]) is gradually deviated from 
its initially Maxwellian distribution. The quasi-linear coefficients of Equation 
(|T5)) are only valid in a small deviation from the Maxwellian VDF of particles 
(Equation ©) and a small damping rate of plasma waves (Equation ([TU)) V Owing 
to these limitations, the analytical method in this paper can only describe the 
initial quasi- linear stage of wave-particle interaction. The ensuing nonlinear stage 
is beyond the scope of this paper, though it can be self-consistently described 
by full particle simulation. 



9. Conclusions and Summary 

Using the linear Vlasov wave theory and quasi-linear resonant wave-particle in- 



teraction for a low-beta multi-species plasma ( Li and Habbal, 2001 Marsch and Tu, 2001[ ),[ 
fast and Alfvcn waves (and their dispersion relations and electromagnetic field 
fluctuations) are analyzed, and then applied to energize a proton/electron/alpha 
plasma with a low proton beta of /3p w 0.01. In this paper, we assume that (1) 
low-frequency Alfven and fast waves are launched from the solar surface with 
the same spectral shape and the same amplitude of PSD, (2) these macro-scale 
waves/turbulences continuously cascade towards kinetic scales at ion cyclotron 
frequencies and ion inertial scales, and (3) the transported kinetic-scale waves 
energize solar wind species through cyclotron and Landau resonances. If the 
wave propagation angle 9 is oblique with respect to the background magnetic 
field, fast-whistler waves can be linearly coupled with hybrid alpha-proton cy- 
clotron waves. Both fast-cyclotron and Alfvcn-cyclotron waves are important 
for preferential heating and acceleration of minor ions. In this paper, noticeable 
energization of alpha particles requires oblique propagation for fast-cyclotron 
waves {0 > 31°) and non-perpendicular propagation for Alfven-cyclotron waves 
{& < 80°). Nearly perpendicular Alfvcn-cyclotron waves at 80° < 9 < 90° signif- 
icantly degenerate their capabilities of energizing minor ions due to very large 
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perpendicular wave numbers at ion cyclotron frequencies. Only through complex 
nonlinear processes, nearly perpendicular Alfvcn waves may heat ions. For in- 
stance, 



Cranmer and Ballegooijen (|2003|) proposed that the damping of kinetic 
Alfven waves successively leads to electron beams, Langmuir turbulence, and 
Debye-scale electron phase space holes, and ions are perpendicularly heated via a 
collision-like scenario. Meanwhile, oblique fast-cyclotron waves have a larger left- 
handedly polarized electric field |(5E+| than the corresponding Alfven-cyclotron 
waves, because the total electric field fluctuation |(5E| of fast waves is stronger 
enough to over-compensate its smaller ratio of left- to right-handedly polarized 
electric fields |(5E+|/|(5E~|. Moreover, in contrast to Alfven-cyclotron waves, 
alpha cyclotron resonance happens at a smaller wave number for fast-cyclotron 
waves especially at « 90°, where the turbulence PSD has a larger local 
amplitude. Hence, fast-cyclotron waves have a greater capability of energizing 
solar wind ions than Alfven-cyclotron waves. Interacting with Alfven-cyclotron 
and fast-cyclotron waves, alpha particles are preferentially energized as the 
primary recipient of wave energy. In contrast to alpha particles, electrons are 
generally less energized. When the differential speed U\\a of alpha particles is 
high enough, Alfven-cyclotron waves lose alpha cyclotron resonance much earlier 
than fast-cyclotron waves. For example, at = 60°, alpha cyclotron resonance 
disappears at U\\a = 0.2 va for Alfven-cyclotron waves and U\\a = 0.6 va for fast- 
cyclotron waves. When alpha particles lose cyclotron resonance, fast-cyclotron 
waves restore to fast-whistler waves without linear mode conversion, and Alfven- 
cyclotron waves begin to mainly energize protons. Moreover, an extremely fast 
beam and/or a high temperature anisotropy of minor ions could release their 
excessive kinetic energy to excite ion cyclotron waves. Therefore both fast and 
Alfven waves are likely to contribute to the eventually sustainable and balanced 
differential speed and temperature anisotropy of minor ions in the solar wind. 

Linear mode conversion between fast- whistler and hybrid cyclotron waves due 
to the presence of minor ions could be an efficient mechanism to dissipate wave 
energy into minor ions. In a warm plasma in the interplanetary space, such mode 
coupling vanishes, and parallel electron heating by Landau damping of oblique 
fast-whistler waves becomes non- negligible. As demonstrated in this paper and 
the previous work of lLi and Habbal1(l200ip . fast-cyclotron waves besides Alfven- 
cyclotron waves could also be considered as a promising candidate for preferential 
heating and acceleration of minor ions in the inner corona. Since fast and Alfven 
waves heat the solar wind ions in very different ways, it is important to know the 
exact nature of plasma turbulence at kinetic scales. In this regard, observation 
data such as Cluster archive can continue to address the issue. 
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Table 1. 


Symbols used 


in this paper and their denotations 


Symbol 




Description 


c 




light speed 


e 




an electron charge 


fcs 




Boltzmann constant 


X, y, z 




Cartesian coordinate system 


Bo, Bo 




uniform background magnetic field along z-axis (-Bq = I Bq ) 


<5B 




magnetic field turbulence 


5E 




electric field turbulence 


j parallel direction (z-axis) 


± 




perpendicular direction (xy-plane) 


P 




polarization of electric field fiuctuation in xy-planc 


+, - 




left- and right-handed electric field polarizations in xy-plane 






unimodular electric field polarization, and 
its left- and right-handed components 


k, fc, Efc 




wave vector (fc = | k|) and its unimodular vector (e^ = k/fc) 


fc||, fej_ 




parallel and perpendicular components of k 


61 




angle between the propagation and the ambient magnetic field 
d = arctan(fcx/fc||) 


i (O. P: 


e) 


plasma species j in an alpha/proton/electron (a/p/e) plasma 


nij {ma 


, rrip, rrto) 


mass of species j 


Zj {Za, : 


^p, 2e) 


electric charge number of species j {za =2, Zp = z^ = 1) 



rij (tIq, rip, ric) abundance of species j 

(neutral charge: ria Za + np Zp + ric Zc = 0) 
(minor heavy ion: n^ <C iip) 





ti^pc) 




^pj ('^P", '^PP) 


plasma frequency LOpj = . Airrij 2? e^/mj 


(He, Op, Qc) 




cyclotron frequency Hj = Zj eBci/{mj c) 
(Qp = 2 He) 


W, Wr, '^i 




complex wave frequency as well as its real and imaginary parts 








«A 


Alfven speed ija = Bo/^iirrip m.p 


(f/||., t/||e) 




differential speed of species j relative to protons along Bp 

(zero current: Ua Za U^\a + l^e Ze l/||e = 0) 


T||j (Tji^, T||p, T||<,) 


parallel temperature of species j 


T±j {T±a, T±p, 


Txc) 


perpendicular temperature of species j 


7j {ic, 7p, 7c) 




temperature anisotropy ■yj = T±j/Ty 



/3j species beta defined as 

(Pa, Pp, /3e) a ratio of species j's thermal pressure to magnetic pressure 

Pj =8 7rnjfcBT||j./Sg 

/3 plasma beta {0 = 0a + Pp + Pc) 
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Table 1. 






""Wi ('^Ik' ""IIP' ''lie) 


parallel thermal speed v^y = ,/2kBTy/mj 






V±j {v±a, V^p, Dxe) 


perpendicular thermal speed v±j = ^2kBTj_j/mj 


fj {/<.,/p, /c) 


gyrotropic velocity distribution function of species j 
in terms of ii||j and v^j 


«i («<.,5p,«c) 


arbitrary-order general resonance factor of species j 


«° «°, ?°, ?°) 


zero-order Landau resonance factor of species j 


«+ «J, «^, «c+) 


first-order cyclotron resonance factor of species j 


"IIj ("11"' "1|P' "lie) 


normalized parallel acceleration rate of species j 


Q\\j (Q||., Qiip, Q||o) 


normalized parallel heating rate of species j 


Q±j {Q±c, Q±p) 


normalized perpendicular heating rate of specie j 



Table 2. Fixed prescription among all initial plasma 
parameters 



Symbol 



Value 



parallel temperature ratio 



4:1:1 



temperature anisotropy of alpha particles 7c 1 



temperature anisotropy of protons 7p 
temperature anisotropy of electrons 70 



Table 3. Parametric study of initial plasma parameters 
Section Figure Propagation Density ratio Differential Proton beta 



angle 1 



speed of alpha 
particles C/iu 



3, 4 


1-8 


0° - 89° 


0.02 







0.01 




5 


— 


0° - 89° 










0.01 




6 


9, 10 


60° 


0.02 


0- 


- 1 ''A 


0.01 




7 


— 


60° 


0.02 







0.001 - 


- 10 
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Figure 1. Dispersion relations of Alfven, fast, and hybrid alpha-proton cyclotron waves at 
{a,b) 9 = 1° and (c,d) d = 89° in a proton/alpha/electron plasma. Note that linear mode 
coupling between hybrid and fast waves happens in panel (c). 
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Alfven Wave 
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Figure 2. For Alfven waves, the dependence of complex frequency {lu = oJr+ii^;) on the paral- 
lel wave number fcii and propagatio 
are indicated by a hatched region. 



lei wave number fcy and propagation angle 9. Heavily damped waves satisfying i.o'i(k) < — ^n 



SOLA: SoPh2011Rev2.tex; 25 February 2013; 23:47; p. 21 



Solar Physics 



a 






a 







0.5 



1.5 



1.0 



0.5 - 







Fas t Wav e 

(a) 




30 



60 



e n 




4 
3 . 

G 

: '^ 3 
: 1 

0.4 
0.3 

0.2 

0.1 
0.0 



a 



90 



Figure 3. For Fast waves, the dependence of complex frequency (ui = uir + ii^O on the parallel 
wave number fcii and propagation angle 9. Linear mode conversion occurs for 6 > 31°, with its 
position in the fcii — 8 domain denoted by a black-and-white line. 
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Figure 4. Properties of quasi-perpendicular fast waves, quasi-perpendicular and quasi-parallel 
Alfven waves in a proton/alpha/electron plasma: (a,e,i) ratio of electric to magnetic field 
turbulence spectra |5Ep/ |5Bp, (b,fj) ratio of left- to right-handed polarizations for electric 
field spectrum |e^| / |e^|, (c,g,k) Landau resonance factors g", and {d,h,l) cyclotron resonance 
factors £+. 
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Alfven Wave 




Figure 5. The dependence of Alfven wave properties in a proton/alpha/electron plasma on 
the parallel wave number fcii and propagation angle 9: (a) ratio of electric to magnetic field 
turbulence spectra |(5Ep/ |(5Bp, (b) ratio of left- to right-handed polarizations for electric field 
spectrum |ej^| / [e^|, (c) Landau resonance factors of protons ^JJ and alpha particles 5^, (d) 
Landau resonance factor of electrons ^^ , (e) cyclotron resonance factor of protons gp , and (f ) 
cyclotron resonance factor of alpha particles ^J . 
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Figure 6. The dependence of fast wave parameters in a proton/alpha/electron plasma on 
the parallel wave number fcii and propagation angle 0: (a) ratio of electric to magnetic field 
turbulence spectra |5Ep/ |5Bp, (b) ratio of loft- to right-handed polarizations for electric field 
spectrum |e]^|/|e^|, (c) Landau resonance factors of protons ^p and alpha particles ^IJ, (d) 
Landau resonance factor of electrons ^J? , (e) cyclotron resonance factor of protons ^p , and (f ) 
cyclotron resonance factor of alpha particles ^ J . 
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Figure 7. The energization capabilities of Alfven waves in a proton/alpha/electron plasma 
on (A) alpha particles and (B) electrons. From top to bottom are (a,d) parallel acceleration an , 
(b,e) parallel heating Qii , and (c) perpendicular heating Qx over the fcii —9 domain. Quantities 
of ail, Qii, and Qx a-fc already normalized by the wave spectrum |(5B(k)p. Here the proton 
energization and the perpendicular electron heating are essentially zero. 
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Figure 8. The energization effects of fast waves in a proton/alpiia/electron plasma on (A) 
alpha particles, (B) protons, (C) electrons. From top to bottom are (a,d,g) parallel acceleration, 
(b,e,h) parallel heating, and (c,f) perpendicular heating over the fcy — 6 domain. Here the 
perpendicular electron heating Q±g is essentially zero. 
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Figure 9. For Alfvcn waves at propagation angle 9 = 60 

f/||^ in a proton/alpha/electron plasma results in different energizations of (A) alpha particles 

(B) protons, and (C) electrons. 
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Figure 10. Depending on the differential speed of alpha particles C/ii^ in a pro- 
ton/alpha/electron plasma, the energization consequence of fast waves at propagation angle 
9 = 60°. Linear mode conversion occurs at f/ii^ < 0.6vai denoted as a black-and-white line in 
the fcii — Uwr, domain. 
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